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We show that graphene monolayer coupled to multi-quantum well system can form a graphene 
spaser, which is a coherent quantum generator of surface plasmons in graphene. The active element 
of graphene spaser is a multi-quantum well system with a design similar to the design of an active 
element of quantum cascade laser. Under realistic parameters of the multi-quantum well system, 
the spasing in graphene monolayer can be achieved at a finite doping of graphene and at a plasmon 
frequency, « 0.15 eV, close to the typical frequency of intersubband transitions in multi-quantum 
well systems. 
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Quantum nanoplasmonics becomes fast growing field 
of nanooptics[l| due to unique property of plasmonic 
system to strongly enhance the optical field within spa- 
tial subwavelength range. Such property results in 
many applications of nanoplasmonic systems, such as 
surface enhanced Raman scattering 0-3, ultrasensingQ, 
surfase-plasmon enhanced photodetectorsjgj], biomedical 
testing^, 0] and others. One of such applications is 
a spaser (Surface Plasmon Amplification by Stimulated 
Emission of Radiation) , which is nanoplasmonic counter- 
part of laser. Spaser has been proposed theoretically in 
Ref. [8] and observed experimentally in Ref. [9]. Spaser 
consists of plasmonic system coupled to active (gain) 
medium, which results in coherent generation of surface 
plasmons with nanolocalized intense optical field. One 
of the limitations on coherent generation of surface plas- 
mons is determined by finite damping rate of plasmon 
excitations. 

In the prese nt p aper we propose to utilize the graphene 
monolayer [Iol - [i"2j as a plasmonic system [l3l - [l6| of spaser. 
Graphene is a layer of crystalline carbon that just one 
atom thick. The honeycomb 2D crystalline structure 
of grap hene results in a very unique energy dispersion 
law[17j. Namely, the low-energy excitations are described 
by chiral relativistic massless Dirac-Weyl equation. The 
chiral nature of electronic states in graphene results in 
strong suppression of electron scattering, which highly 
increases the mobility (~ 15,000 cm 2 V _1 s _1 ) of charge 
carrier in graphene and strongly suppresses theoretically 
estimated plasmonic losses. For plasmon frequency be- 
low the interband loss threshold, which is determined by 
the level of graphene doping, i.e., the Fermi energy, and 
below the optical phonon loss threshold, w 0.2 eV, the 
plasmon relaxation time can be estimated from DC mea- 
surements and is around r rs 0.6 ps [l4j |. which results in 
low plasmon damping rate 1/ujt rs 0.01. 

As an active element of graphene spaser, we will con- 
sider a multi-quantum well system. The design of such 
system is similar to the design of an active element of 
quantum cascade laser. Quantum cascade lasers, pro- 
posed theoretically in 1971 [18j and realized experimen- 
tally in 1994,[lj|[20| utilize the intersubband optical tran- 
sitions in specially designed superlattice of quantum wells 
to generate coherent lasing in mid-infrared and terahertz 
optical range. [ll| Within the active region of quantum 



cascade structure, optical transitions occur between sub- 
band levels of dimensional quantization in the growth 
direction. [22j The advantage of quantum cascade laser 
system over conventional semiconductor lasers, which are 
based on inter-band transitions, is that a singe electron 
propagating through the cascade structure can generate 
as many photons as the number of active regions within 
quantum cascade laser. In the graphene spaser we utilize 
only one multi-quantum well active element of quantum 
cascade laser. The design of such active element is char- 
acterized by strong electrical injection of electrons into 
the upper subband level of multi-quantum well system 
and fast and efficient extraction of electrons from the 
lower subband level. 

Our spaser system is shown schematically in Fig. [1] 
for one of the possible designs of active element of cas- 
cade system. Since the plasmonic modes of graphene 
are strongly localized in the growth direction, i.e. along 
axis z, with localization length « 3 — 5 nm, only one 
active element of quantum cascade system is coupled to 
the plasmonic modes of graphene. The gain, i.e., popu- 
lation inversion, in multi-quantum well system is intro- 
duced through a voltage applied to the system. 

To find the plasmonic modes of graphene, we assume 
that graphene monolayer is surrounded by two homoge- 
neous effective media with dielectric constants ej and t^. 
Therefore we replace the cascade system with homoge- 
neous dielectric media. Spatial distribution of electric 
field of plasmonic excitations in graphene and their dis- 
persion relation are determined by solution of Maxwell's 
equations with corresponding boundary condition. The 
boundary conditions are introduced in terms of 2D cur- 
rent, generated in graphene layer, where the current is 
proportional to the in-plane electric field and is deter- 
mined by graphene conductance, a{ui,q). We consider 
TM plasmon mode in graphene monolayer with the dis- 
tribution of electric field in the following form[l4j] 
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FIG. 1: (a) Schematic illustration of coupled graphene and 
quantum cascade laser structure. The numbers (in nm) are 
the well (Gao.47Ino.53As) and barrier (Alo.4sIno.52As) widths. 
The wells and barriers are shown by light and dark gray col- 
ors, respectively, (b) Energy-band diagram (schematic) of the 
active region of a quantum cascade laser structure under con- 
stant electric field. The relative wave functions of subbands 
(1) and (2) as well as optical transitions corresponding to 
spaser action are shown schematically. 



at z > 0. The corresponding profile of magnetic field is 
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Here £ q is the amplitude of the plasmonic mode with 
frequency oj q and wave vector q. The spatial decay of 
electromagnetic field away from the graphene layer is de- 
termined by wave vectors k± g and k 2z , which are given 
by the following expressions 
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Considering the limit q 3> lo/c and assuming the Drudc 
expression for graphene conductance, we obtain from the 
boundary conditions the well-know dispersion relation for 
plasmons in graphene 



q(u) « k lz « k 
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Ae^Ep v uu q , 

where Ep is the Fermi energy of electrons in graphene 
and r is the relaxation time, which also determines the 
plasmon decay time. 



Applying the second quantization approach to plas- 
mons in graphene, we introduce electric field operator in 
the following form 



E 
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where a q and a~j~ are operators of creation and annihila- 
tion of a plasmon with the wave vector q and frequency 
uj q . Within the second quantization approach, the ampli- 
tude E q of plasmonic mode is determined from the con- 
dition 
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where So is the area of graphene layer. 

The Hamiltonian of a spaser, shown in Fig. [TJ has the 
form 



hJ2n 12 {k,q) (a+c+ k c 2 ,k + h.c.J , (10) 

q ,h 

where the first sum describes the Hamiltonian of the gain 
medium, which consists of two subbands "1" and "2" 
(see Fig. [TJ with energies E\(k) = e\ + h 2 k 2 /2m* and 
E 2 (k) = e 2 + h 2 k 2 /2m*, respectively. Here m* is an elec- 
tron effective mass, e\ < e 2 , and k is a two dimensional 
electron wave vector. The second sum in Eq. (fT0|) is the 
Hamiltonian of plasmons in graphene layer, and the last 
sum describes the dipole interaction of plasmons with 

active medium, cx Ed. Here the dipole vector d of the 
active medium has only z component. The dipole inter- 
action strength is determined by the corresponding Rabi 
frequency, fii 2 (g), 



n 12 (k,q) = e £ q / drip^ k (r)E qiZ zip 2tk (r), 
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where eo is the electron charge, ipi t k and ip 2 ,k are elec- 
tron wave functions with wave vector k of subbands 1 
and 2, respectively. Assuming that the electron dynam- 
ics in x, y, and z directions is saparable, we can ex- 
press the wave functions in the following form i/^.fc = 

— 1/2 "* - 

Xi(z)S exp(ikp), where p = (x,y) and i = 1,2. Then 
expression (fTTj) becomes 
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where d\ 2 = J dzx*{z)eozx2(z) is the dipole matrix ele- 
ment, and L is the distance between graphene monolayer 
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and the maxima of the subband wave functions Xi( z )i 
which correspond to the active quantum well of cascade 
structure. 

Introducing the density matrix of the active region of 
the system, p%{t), we describe the dynamics of the spaser 
within the density matrix approach, writing the equation 
for the density matrix, ihpAt) = [Pj:{t),'H], and for oc- 
cupation numbers of the plasmons. Within a stationary 
regime we obtain the condition of spasing in the following 
form! 



(7 + ri 2 ) 5 
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>7ri2, (13) 



where 7 = H/t is the surface plasmon relaxation rate, Ti2 
is the subband polarization relaxation rate, and o>i 2 = 
(e 2 —£x)/fl. Substituting expression ([9]) into Eq. (fT2|) . we 
obtain 
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Here the wave vector fc , i.e., quasi- Fermi wave vector, 
determines the largest wave vector of the states, which 
are partially occupied by electrons in the active region. 
Only such states contribute to the sum in Eq. [13] 

At the resonance, i.e. at u\ 2 = u qi the condition of 
spasing becomes 
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For a given structure of the active region, expression f| 15[) 
introduces a constrain on plasmon relaxation rate, 7. 
The typical parameters of the active region (correspond- 
ing, for example, to the structure shown in Fig. [T|) arc 
hcui2 = 130 meV, d\ 2 = 1.5eo (nm-eo), and L ~ 16 nm. 
The typical intersubband relaxation time in th e q uan- 
tum well cascade systems is t 12 = fi./r 12 ~ 1 ps[23l l24j. 
Under resonant condition the inter-subband transition 
frequency Huji 2 = 130 meV, which is equal to the plas- 
mon frequency, correspond to the plasmon wave vector 
q^ 1 m 30 A- see Eq. ([5]). Assuming that k ss 0.04 nm -1 
and ei + e 2 = 5, we obtain 



(16) 
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7 < 2.5 • 10" J eV 

with corresponding relaxation time r > 3 ■ 10 

Substituting relation © into inequality (fTS"]) , we obtain 
a condition on plasmon relaxation time 
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The dependence of the minimum plasmon relaxation 
time, r m j n , on the plasmon frequency w q is shown in 
hg- Ula) for different Fermi energies Ep of electrons in 
graphene. The data clearly show that there is a broad 
range of plasmon frequencies, hw q ps 0.1 -j- 0.3 eV, at 



which the minimum relaxation time is relatively small, 
Tmin ~ 0.03 ps. Outside this range the minimum re- 
laxation time strongly increases. The typical values of 
frequencies of intersubband transition in quantum well 
structures are within the frequency range 0.1 -j- 0.3 eV. 

The condition (| 1 T[) can be also expressed in terms of 
dimensionless plasmon damping rate, 7 = l/(cj g r), 



7 < 7n 
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Restriction (|18[) on the surface plasmon damping rate 
depends in the plasmon frequency, uj q , which under res- 
onant condition is equal to the frequency of intersub- 
band optical transitions within the active region of cas- 
cade system. As a function of frequency u) g , the damp- 
ing rate 7 ma x has a sharp maximum, the position of 
which is determined by condition q « 1/2L or uj q w 

1 /2 

(2e 2 Ep /[LH 2 (ei + e 2 )]) , and depends on parameters 
of the system, such as the Fermi energy, Ep , of electrons 
in graphene and the distance L between the graphene 
monolayer and the active quantum well of cascade struc- 
ture. The dependence of the plasmon damping rate 7 ma x 
on the frequency of surface plasmons, ui q , is shown in Fig. 
[2jb) for different values of the Fermi energy Ep. The 
maximum value of 7 max is fc 2 |di 2 | 2 /[4eFi 2 (ei + e 2 )] « 0.2 
and does not depend on the Fermi energy. The Fermi 
energy affects only the position of the maximum. 

Below the frequency of optical phonons in graphene 
monolayer, which is ~ 0.2 eV, and below the Fermi fre- 
quency, which determines interband loss threshold, the 
relaxation time of surface plasmons in graphene can be 
estimated from the DC conductance of graphene. Such 
estimations show that the plasmon relaxation time in 
graphene is ~ 0.6 ps, which corresponds to the plas- 
mon damping rate 7 « 0.01 at frequency Hu) q = 130 
meV. Such plasmon damping rate is well below 7 max for 
a broad range of frequencies uj q - see Fig. . 

In Ref. [l6j the plasmon damping rate 7 = l/(ui q T) in 
graphene has been measured experimentally by infrared 
nano-imaging method. The best fitting to experimental 
data has been obtained for strong damping 7 = 0.135 
at plasma frequency hijj q — 0.11 eV, which corresponds 
to the plasmon relaxation time r = 0.044 ps. Even at 
such strong plasmon damping, the plasmon damping rate 
7 = 0.135 ps is less than the maximum rate 7 max ~ 
0.2 - see Fig. HJb). Taking into account that plasmon 
relaxation rate strongly increases above optical phonon 
loss threshold, which is f» 0.2 eV, we can conclude that 
the optimal Fermi energy in graphene should be around 
0.5 eV. For such Fermi energy the maximum of 7 max is 
realized at plasmon frequency huj q ks 0.15 eV - see Fig. 
[5Jb) . The above results clearly suggests that the spasing 
can be realized in realistic graphene systems coupled to 
quantum well structures. 

In conclusion, the graphene monolayer deposited on a 
top of specially designed stack of quantum wells forms 
a spaser, which is a coherent source of surface plasmons 
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in graphene. The system of quantum wells is designed 
similar to an active element of quantum cascade laser 
and is characterized by strong electrical injection of elec- 
trons intro the upper subband level of multi-quantum 
well system and fast and efficient extraction of electrons 
from the lower subband level. Under realistic param- 
eters of the multiple quantum wells, our analysis show 
that the condition of spasing, i.e., coherent generation 
of plasmons in graphene, can be achieved at finite dop- 
ing of graphene monolayer, which corresponds to elec- 
tron Fermi energy in graphene around 0.5 eV. The above 
analysis can be applied not only for graphene monolayer, 
but also for graphene-based nanoscale structures such as 
graphene nanoribbons and graphene quantum dots. 



FIG. 2: (a) Minimum relaxation time, defined by Eq. (|17|l . 
of plasmons in graphene is shown as a function of resonant 
plasmon frequency for different Fermi energies of electrons 
in graphene. The spasing is realized at plasmon relaxation 
time t > T m in. (b) Plasmon damping rate 7 ma x, defined by 
Eq. (|18|l . is shown as a function of the plasmon frequency 
ui q . The dashed line shows the experimental value of plasmon 
damping rate 7 « 0.135 [lg|. The numbers next to the lines 
in (a) and (b) are the Fermi energies in eV. 
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